IP/BBSR/2003-01 



Nonsingular Cosmologies from Branes 



Anindya Biswas, Sudipta Mukherji and Shesansu Sekhar Pal 

Institute of Physics, Bhubaneswar- 751 005, India 

anindyab, mukherji, shesansuOiopb .res . in, 



We analyse possible cosmological scenarios on a brane where the brane acts as a dy- 
namical boundary of various black holes with anti-de Sitter or de Sitter asymptotic. 
In many cases, the brane is found to describe completely non-singular universe. In 
some cases, quantum gravity era of the brane-universe can also be avoided by prop- 
erly tuning bulk parameters. We further discuss the creation of a brane-universe by 
studying its wave function. This is done by employing Wheeler-De Witt equation in 
the mini superspace formalism. 
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1 Introduction 



Motivated by string theory, the AdS/CFT correspondence and the hierarchy problem in particle 
physics, brane-world models have been a focus of interest in recent years - (Hj- In these models, 
our universe is realized as a boundary of a higher-dimensional space-time. In this context, a well 
studied example is when the bulk is an AdS space. The gravitational interaction among matter 
on this brane is found to be described by standard laws when one considers distance scale much 
larger than the AdS length scale 

In the cosmological context, many authors have considered (see for example |5] - |28p 
the embedding of a four dimensional Friedmann-Robertson- Walker (FRW) universe in five di- 
mensional bulk geometry. The bulk is often described by AdS or AdS Schwarzschild black hole 
metric. In the later case, the mass of the black hole is found to act effectively as an invisible 
energy density on the brane with the same equation of state of radiation matter. This, however, 
is found to have a nice interpretation in terms of AdS/CFT correspondence. The correspondence 
allows one to reinterpret the AdS-Schwarzschild geometry as a source for a four dimensional 
CFT at finite temperature. As pointed out in - it is easy then to understand why the 
radiation-dominated FRW universe emerges. This is because all CFTs have the same equations of 
state up to numerical constants and hence the FRW equation takes the same form as in the case 
of radiation dominated universe. Replacing the bulk now by a five dimensional AdS-Reissner- 
Nordstorm black hole may seem like a straight forward generalization. However, as it was found 
in |32j . such geometry has far reaching consequences on the brane evolution. First, FRW equation 
now describes the brane as dominated by induced radiation and "stiff" matter. However, in the 
equation, the stiff matter contribution comes with an opposite sign*. As a result, the universe 
(flat, open and closed) is never found to shrink to zero size. It is important to note that the bulk 

*This was first noticed in |13|. 
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electric charge plays the role of a "regulator" on the brane. As we take the charge to zero, big 
bang singularities are found to reappear^. 

Here, in this paper, we further analyse the D dimensional brane evolution in various D + 1 
dimensional AdS or dS black hole backgrounds of various topologies. We first study these back- 
grounds carefully and find critical values of mass and other parameters in order these geometries 
to correspond to black holes with non-degenerate horizons. We believe that some of these explicit 
results here are new. We then study the brane dynamics in these geometries. It is found to be 
described by FRW equation which has contributions from induced radiation matter, stiff matter 
and an effective cosmological constant. This constant depends on a certain combination of bulk 
cosmological constant and brane tension. We then analyse the classical dynamics of the brane by 
solving the FRW equation in case of vanishing and nonvanishing effective cosmological constant. 
We find that all open, flat and closed universes have no big bang singularities when the bulk 
has a nonvanishing electric charge. For the closed universe, the size is bounded from above and 
from below. On the other hand, for open and flat universe, the brane starts with finite radius 
and expands for ever. We also find here that the minimum radius of the universe can be made 
sufficiently large by tuning the bulk mass and charge to large values. This, in turn, means that for 
such universes we can avoid reaching a quantum gravity era in the past. This scenario, in a more 
general set up, is analysed recently in |34| . We then analyse the brane universe when there is an 
effective (negative or positive) cosmological constant on the brane. Though the FRW equations 
here can not be solved analytically in a closed form, we study the equation by numerical means. 
The non-singular nature of the solutions is also found to persist in this case. However, is all the 
above cases, universe is found to be singular when the charge is set to zero. We, therefore, turn 
our attention to the fate of these singularities from the perspective of "quantum cosmology" . By 
employing the Wheeler De Witt (WDW) equation in the mini superspace formalism, we find that 
the minimum size of the brane universe is of the order of the Planck length. We feel that this 
may indicate that the universe stabilizes against collapse due to quantum effects. However, we 
may add here that quantum cosmology at this present stage has certain inherent ambiguities. 
These are associated with the orderings of operators, choice of boundary conditions of the wave 
functions etc. That is why our results here will be more of speculative in nature. 

This paper is organized as follows. In the next section, we critically analyse various topo- 
logical AdS and dS black hole backgrounds. Considering brane as a dynamical boundary of these 
geometries, we set up brane equation of motion by studying the brane effective action. In section 
3, we analyses classical evolution of the brane by explicitly solving brane equation of motion, $ 
We then critically analyse various physical properties of our solutions with a particular emphasize 
on their early and late time behavior. Section 4 of this paper deals with the brane dynamics from 
the WDW perspective. We solve relevant WDW equations , in mini superspace, for tensionless 
branes and analyse its behavior at short distance scales. For branes with tension, one needs to re- 
sort to some approximate methods. We study these branes in section 5 via WKB approximation. 
The paper ends with a discussion of our results. 

^Sortie earlier work on brane cosmology can be found in 32 33 . 
"'"Some cosmological scenarios on the brane have been reviewed in |27|. 
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2 Brane effective action 



Bulk with AdS asymptotics 

The dynamics of a D dimensional brane inside a bulk gravitational field is dictated by the 
following bulk-boundary action with a negative cosmological constant A = — gjs ^ 

1 f ^+1 /—( r> „2 , D(D-l) 



f d D+1 xMR-F 2 + 
Jm v 



where 



IQttG d+ i Jm \ I 2 

87TGd + i JdM 07TGd+1 JdM 

where F 2 is the Maxwell kinetic energy term and K, is the trace of the extrinsic curvature K,^ 
taken with respect to the induced metric on the brane 7^ = g^ u — n^n u . This term is required for 
a well-defined variational principle on the space time boundary [SSI- Here is the unit normal 
vector to the brane. Furthermore, T in Q is the brane tension and this term is required in order 
to get a finite action and stress tensor jSHlEZl^- Variation of (0) with respect to the induced 
metric gives us 

T 

K-viv = D _ ^ iw ( 2 ) 

A large class of solutions of bulk equations of motion can be represented by 

da 2 

ds 2 = -h(a)dt 2 + -j-j-j + a 2 jijdx l dx j , (3) 

a 2 lo d M {D-1)lo 2 d Q 2 16vrG D+ i m 

h{a) = k+ p-^ + 8(D-2)a 2 ^ WD= (2?-l)Vb-i- (4) 

Here, k = 0, =Fl, correspond to flat, hyperbolic and spherical geometries of D dimensional subspace 
for a given a. ^ij is the metric for a constant curvature manifold M B_1 with volume Vd—i = 
Jci' D_1 Xv / 7- As can easily be seen from and (jU, for M = Q = 0, the bulk is a simple D 
dimensional AdS space. For only Q = 0, the bulk is an AdS Schwarzschild black hole, while for all 
the parameters non-zero and within certain domain, the background corresponds to charged AdS 
black holes. The parameters M and Q can then be identified with Arnowitt-Deser-Misner mass 
and charge respectively. Some analysis of the causal structure of this metric has been performed 

in maun!- 

We would like to comment on the horizon structure of the metric (jHJ). First, note that this 
metric asymptotically is anti-de Sitter for all values of parameters. However, for the metric to 
describe the exterior of black hole with non-degenerate horizon, we need to restrict the parameters. 
To have such a geometry, we would like a^ 2D ~^h(a) to have a simple root at a = <iq such that 
h{a) is positive for a > ao- This gives a lower bound on the parameter M. This can be seen as 
follows. Let h(a) be zero for a = ao. Therefore, we can invert the relation to get mass M as a 

§In general, one should add more derivative terms to see |37| . 
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function of ag and other parameters as 



M = ^ 



(D - l)u 2 D Q 2 a 



8(D-2) 



(5) 



Now, by taking first and second derivative of © with respect to ao, we see that M has a minimum 
M crit at og rit when 

8D(a§ rit ) 2{D_1) + 8(D - 2)kl 2 {af t f D ~ A - l 2 u 2 D (D - l)Q 2 = 0. (6) 

M cnt can therefore be determined by first solving ag from @ and then substituting it into 
(JsJ). If we now define a#(M, Q) as the larger of the positive solutions, we will have the following 
inequality: 

af\Q) < a H (M, Q) < oo, as M crit (Q) < M < oo. (7) 

Though, in higher dimension, expressions of M cnt and <2g nt are either not illuminating or hard to 
obtain, but for D = 3, we have 



M 



crit 



crit 

a 



V2(-k 2 l 2 + 3tu 2 Q 2 + fcZ^/fc 2 ! 2 + 3u;|Q 2 



\ 



(8) 



We note here that aQ nt — * as Q — > 0. For D = 4, that is when the bulk is five dimensional, it is 
also possible to find M crit and aQ rit explicitly. However, the expressions are very big except for 
k = 0. In this case, 



Af 



crit 



33 (huoiQ 2 + y^wfQ 4 
16(/ 2 a; 2 Q 2 + ^/VQ 4 
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crit 

a 



1 

36 



y(zVg 2 + v //MQ 4 )«- 



(9) 



As we would like to find the dynamics of a D dimensional brane moving in (|3"]). it is instructive 
to find out the effective action that controls the dynamics. For this purpose (and hence in the 
rest of the section), we do not need to explicitly specify the functional form of h(a). As we will 
see, effective action describing the boundary dynamics can be expressed in terms of h(a). 

The extrinsic curvature, /C^„, can be calculated following Those are 



and 



/CI 



y/h(a) + a 2 



h (a) + 2d 
2^h(a)+a 2 ' 



(10) 
(11) 
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where h (a) 



dh(a) 
da 



and % = 1, . . . , _D — 1, #j's are the angular coordinates, a 



^ and r denotes 



the proper time as measured along the brane world volume. Details of these computations can be 
found in many earlier literature, see for example jlU]. Inserting these values of extrinsic curvature 
in (J2J), we get 

= —v (12) 



y/h(a) + a 2 



The dynamics of the brane can be captured through an effective Lagrangian L. Clearly, 
the only degree of freedom in L will be a - so the Lagrangian which we are about to write down 
is a construction in the mini-superspace. It is given by 



L = M, n 



a sinh 1 



VHa) 



h(a) + a 2 + 



D - 1 



(13) 



where M p is Planck mass, a dimension full constant. The Euler-Lagrange equation for a then is 
given by 

h'(a) + 2a _J_ 
2^h{a)+a 2 D-l 1 ' 

L then reproduces the correct equation of motion of the brane @ along with ()12|) . We see from 
the Lagrangian, that the momentum conjugate to a is 

p = ^ = M p sinh- 1 . (15) 

da F [vH a ) 

It is then straightforward to construct the hamiltonian, 

T 



H = pa - L = M p 



h(a) cosh f — — 



D-l 



(16) 



We will, in section 4, use this hamiltonian with suitable ordering to analyse a possible quantum 
scenario for the brane universe. However, before doing so, we will study the classical equation of 
motion of a in the next section. 

Bulk with dS asymptotics 

Brane evolution in de Sitter black hole background can simply be obtained from (j3J) by 



replacing V 



-I and setting k = 1. In this case, h(a) = can have three positive roots. The 



largest one corresponds to the cosmological horizon and the rest two correspond to the inner and 
outer horizons of the black hole |41| . The brane equations remain same as in equations ()12|) and 
(|14|) with the hamiltonian as in ()16|) . 



3 Classical cosmological evolution of the brane 



The classical cosmological evolution of the brane is determined from the solutions of Q12j) . Sub- 
stituting the form of h(a) from Q in (|12|) . we get 

k u D M (D - l)u;lQ 2 

a* + A ° + ^ " 8(D - 2)a^~2 » ( 17 ) 
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where Ad = rg^zrns ~~ jz is the effective cosmological constant on the brane. Note that for AdS 
bulk, the effective cosmological constant can be either zero, negative or positive. However, for dS 
bulk, Ap is strictly positive. 

We can now find solutions of (|17|) in different situations, namely, for different values to 
k, for different choices of Ad (positive, negative and zero). We thus have three distinct cases: 
(1) Q = 0, M = 0, (2) Q = 0, M > and (3) Q > 0, M > and for each of these three cases one 
have the following 9 situations 

A D = 0, fc = 0,±l 
A D >0, fc = 0,±l 
A D <0, fc = 0,±l. 

(18) 

We will now discuss all these cases in the following. 
(1) M = Q = 

For M = Q = 0, the brane universe moves in AdS/dS bulk. Solving (|17|). we get the time 
dependence of the scale factor of the brane as 

1 -.e /KEr for k = 0, 



A D 

=— coshf VAdt ) for k = +1, 



= !— sinh(\/Aor) for = —1. 
VAd v y 

(19) 

The above solutions are valid only for Ad > 0. For k = 0, we get a de-Sitter like expansion. 
For k = 1, the universe initially is infinitely large, then it shrinks at finite time to a finite size 
controlled by \/Ad and, subsequently, expands and reaches infinite size at late time r. For 
k = — 1, the radius starts out from zero size and then grows to infinite size with time. 

For the case of Ad < 0, there are no physically relevant solutions except for k = — 1, where a 
behaves as 



-^sin(V|A DW . (20) 



(2) M > 0,Q = 

In this case, the brane moves in Schwarzschild black hole with AdS/dS asymptotics. The effective 
brane cosmological constant A^ can be either zero or non-zero. We discuss both these cases for 
the D + 1 bulk in the following. Though in many cases exact solutions can be obtained, we resort 
to numerical analysis where exact results are either hard to obtain or less illuminating. 

(i)A D = 
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It turns out that, when the effective brane cosmological constant is zero, the exact solutions are 
most easily obtained in terms of conformal time 77 defined as a(j])dr] = dr. The solutions of (|17jl 
then gives, 



0(77) 



D 



■\jMoj D rj 



2 

D-2 



for k = 0, 



1 D-2 2 
(MLU D )D-2 S - m [ — — 77] d-2 for A; = +1, 

1 £) - 2 2 
(M(j£))D-2 S i n h[ 77] -°- 2 for & = -1. 



(21) 



As can be seen from above, for k = 0, the brane starts out from singularity and expands to 

infinite size at late time. However, the rate of expansion and acceleration depend crucially on the 

dimension of the brane. For D = 4, it expands at a constant rate and for D > 4, the expansion 

rate decrease with time 77 producing a deflationary scenario. In terms of cosmic time r, the 

brane expands as radiation dominated universe a ~ td. For k = 1, we get cyclic universe with 

1 

maximum size of the brane is dmax — (Mujd) d - 2 where M is the mass of the bulk black hole. 
For k = — 1, on the other hand, we have an open universe as expected. 

(ii)A D + 



When the brane has effective cosmological constant (negative or positive), the dynamics of the 
universe can be analysed in a quite simple manner. We discuss first the case of positive cosmo- 
logical constant and then the case with < 0. 



{a) A D > 



In this case, (|17jl can be solved exactly in arbitrary D for k = with the result, 



(r) = ^— ^— ^ D sinh[ — — r] o for k = 0. 



(22) 



Therefore, the brane here emerges from the bulk singularity and then expands for ever. In 
particular, at late time, we have de-Sitter like inflation as can be seen from expanding (|22j) for 
large r. 

However, for k = ±1, for only D = 4, exact analytical result can be found and is discussed for 
k = 1 in some detail in |42| I43j . For other dimensions we could get the behavior of the scale 
factor through numerical analysis. For D = 4, we get 



a(r) 



2A 4 



for 4A4W4M = 1, 



V2A 4 



V 7 ! - 4cj 4 MA 4 cosh(2VA~4r) + k] 2 for 4w 4 MA 4 < 1, 
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y/2A 4 



j. 

v/iu^MAl^T sinh(2v / A 4 r) + k 2 for 4cj 4 MA 4 > 1. 



(23) 



Thus, for the case 4A4U4M = 1, we get for de-Sitter expansion for all values of k for large r. 
However, for small r, expansion is sensitive to k. For k = 1 the brane starts from a finite size 
, while for = — 1, brane again emerges from bulk singularity. As can be seen from (|23[). 

when 4A4W4M < 1, the universe radius is bounded from below by — 4A4UJ4M + k)2 /^2A 4 for 
fe = 1. 



As mentioned earlier, for D > 4, analytical solutions in closed forms are hard to obtain. However, 
below are the plots of a with time for different k. 



a 




31 . 63 [ , , I 

0.02 0.2 0.4 

n 

Figure 1: The universe for D = 6 and k = 1, loqM = 100, Aq = .001. It begins at a finite size 
and expands for ever. Here and in other figures, all the quantities are measured in appropriate 
powers of length. 
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Figure 2: The universe for D = 6 and uqM = 100, Ag = 1. Here k = —1 corresponds to the 
solid line while dotted line is for k = 0. In both cases, the universe begins from singularity and 
expands for ever. 



(b)A D < 



The case when effective cosmological constant is negative on the brane, (|17|) can be exactly 
integrated for arbitrary D for flat universe with the result 

/Mu D \jj . .DyfkE 2. 

a ( T ) = ( ~~ a — ) sm i T \ * or ^ = 0- 

\ Ad ' 2 

(24) 

We therefore have cyclic universe with minimum radius being zero while the maximum is deter- 
mined by the ratio of black hole mass and cosmological constant. 

As before, for k = ±1, exact solution can be found easily for D = 4. Those are given by 

i 

2 for k = ±1. (25) 

For D > 4, the behavior of a can be found by numerical means and is shown in the following 
figure. 

(3)M > 0, Q > 



a(r) 



y/AMu 4 A 4 + 1 sin(2v / A7r) - k 



In this subsection, we turn our attention to the case where the brane is moving in charged AdS 
background parametrized by mass M and charge Q. In this case, as will be discussed below, we 
get a host of brane cosmologies which are completely non-singular during the brane evolution. 
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L 

0.4 0.8 1.2 1.6 

z 

Figure 3: Evolution of the brane in D = 5 for u$M = 500, A5 = .01. The dotted, dashed and 
solid lines correspond to k = 0, 1, —1 respectively. 



(i)A D = 



We start here with the case, where the effective cosmological constant is zero. In this case, 
Q17|) can be exactly solved for any D with the result 



a{rj) 



1 AD-l)Q 2 u 2 D . M 2 u 2 
Mu D V 8{D - 2) 



+ 



D 



(D - 2)W 



Muj d 



2 

Mlo d 



1 + 



2(D - 2) M 2 



sin((L> - 2)r? 



D 2 for k = 0, 

for = +1, 



1 

D-2 



1 + V 1+ 2^)S' :OSh « D - 2) '' 



1 

D-2 



for 



(26) 



It is important to note that k = 1 solution only makes sense when the charge satisfies the 
inequality 



Q < v / 2M 1 



I D-2 
D-l 



(27) 



For D = 4, these configurations are studied in some detail in [22]. As can be seen from the figure 
4, all three solutions are non-singular. For k = 0,-1, universe starts from nonzero value of a and 
expands forever. Thus for flat brane 



1 (D - l)Q 2 m 2 D 
Mu D 8(D - 2) 



2, .2 . 1 

D-2 



< a < 00, 



and for k 



Mu D 



1 + 1/1 + 



(g - i)Q 2 

2(L> - 2)M 2 



1 

D-2 



< a < 00. 



(28) 



(29) 
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On the other hand, we have bouncing non-singular universe for k = 1 where the universe oscillates 
between two non-zero values of a as 



{—\}-\l 1 -2(D-2)M*)\ - fl -lTrf"2(Fp)l • (30) 

Notice that as we take Q to zero, the minimum radius collapses to a singularity. So, in a sense, 
the charge parameter Q acts as a regulator. This behavior is also intuitively expected. As in this 
case, the brane moves in electrically charged background, the flux lines have to end on the brane. 
These fluxes, in turn, do not allow the brane to shrink to zero size. A typical behavior of the 
universe for k = 0, ±1 is shown in the following figure. 



14 - 

a 

8 - 

2 - 
- 




1 2 3 4 5 

1 



Figure 4: The universe for k = —1,0, 1 are shown in solid, dotted and dashed lines respectively 
for D = 5, oj^M = 500, 0J5Q = 1 and A5 = 0. As can be seen, in all these cases, big bang like 
singularities are absent. 
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We further would like to emphasize that for k = — 1 and k = 1, the minimum size of the 
universe can be made sufficiently large by taking bulk mass and charge to large values. Thus, 
one would expect to have quantum gravity corrections to be small even when the brane is at its 
minimum size. In this sense, the universe avoids the usual quantum gravity era in the past. 

(ii)A D ^ 

It turns out that in this case, it is hard to integrate (|17j) to get a closed form expression for 
a as a function of r. However, the equation can easily be integrated numerically. In what follows, 
we discuss the solutions for Ap > as well as < for D = 4. 

{a)A D > 

As can easily be checked from (|17|) that for fixed Ap, for certain range of parameters 
(M, Q), a becomes zero at a positive value of a at some finite time r for all k. The universe then 
begins at that radius. At late time, brane expands exponentially independent of the values of k. 
The expansion rate of course depends on the effective cosmological constant Ap. The complete 
evolution of the brane is shown in the figure below for certain choice of parameters. 




Figure 5: This shows typical behavior of the universe for A4 = 1. Here again k = —1,0,1 are 
shown in solid, dotted and dashed lines for D = 4, W4M = 50, oj^Q = 0.1. Universe begins at a 
finite size and expands at late time. 
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(b)A D < 



In the case where the effective cosmological constant is negative, for fixed A^, for certain 
choice of parameters (M, Q), there are two real positive roots of the equation a = for all k. The 
universe then begins then at finite a m { n and ends at another value of a max . Below a m i n and above 
a max , the brane is unstable. The complete evolution of the brane for certain choice of parameters 
is shown in the figure below. 



3 

2 . 5 

2 

a 1.5 
1 

. 5 

-0.4 -0.3 -0.2 -0.1 

?7 

Figure 6: This shows typical behavior of the universe for A4 = — 1. k = — 1,0, 1 are shown in 
solid, dotted and dashed lines for D = 4, 004M = 100, L04Q = 0.1. Universe begins at non-zero 
size and ends at finite size. Instabilities arise below the minimum and above the maximum radii. 

As we have now seen, in some instances, universe begins at the singularity of the bulk. This 
in turn imply big bang like singularity on the brane. In the next section, we study these singular 
universes from the perspective of "quantum cosmology" . 



4 Quantum cosmologies arising from Wheeler-De Witt Equation 

In this section, we turn our attention to a quantum version of the cosmological models of 
previous section. Due to lack of an existing formulation of quantum gravity, our task becomes 
really difficult. We, therefore, need to make drastic simplifications. The approximation that is 
most commonly used is the mini superspace truncation of the Wheeler-De Witt formalism. In 
what follows, we will make such approximation in order to discuss a suggestive quantum version 
of our previous models. Even with in the mini-superspace approximation, quantum cosmology 
suffers from many ambiguities. The operator ordering ambiguity in the hamiltonian operator, 
choice of boundary condition on the wave function are among many others. We may say that 
we have nothing to add to improve our understanding on those issues. We rather work with 
a particular operator ordering and chose a particular boundary condition on the wave function 
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which seems natural to us. As a result, our result in this section would at most be speculative in 
nature. We will closely follow the method developed in 01]. Creation of brane-universe in AdS 
space is discussed in [Tl*] . 



-*35 and is 



First, the Wheeler-De Witt hamiltonian can be obtained from ()16|) by substituting p — > 

(31) 



H = M p 



(J_d_\ ,1 T 



h* cos — — h 4 

\M p da) (D-l) 



and the evolution of the brane universe is then governed by 

Hip{a) = 0. (32) 

Here, ip{a) is the wave function of the brane universe while H is a quantum hamiltonian of the 
brane. In writing down (|31|) . one faces the ordering ambiguity. Hermiticity of the hamiltonian is 
used as guideline to set this operator ordering. We, however, may note that the ordering here is 
not unique and one can, in principle, work with a different ordering prescription. This, in turn, 
will change the nature of the quantum behavior. 

It is now straightforward to solve 1)32(1 particularly for T = 0. The solution is 

i>mn{a) = C mn ((j) m - 4> n ), (33) 

with 



i 



4>j = h * exp 



(j+~)xaM p 



(34) 



C mn here are normalisation constants which can be explicitly evaluated. However, for our purpose, 
we would not require their explicit forms. In the above expression of the wave functions, m, n 
are the integer valued quantum numbers describing the state of the brane. Negatives values of 
m, n do not give normalizable wave function so we discard them. Same is the case for m = n. 
It is now quite easy to see that for both M = and Q = 0, and M^0,Q/0, the solution, as 
written above, satisfies both the boundary condition (given later) and Wheeler-De Witt equation. 
However, for M ^ 0, Q ^ 0, we have another solution, which also satisfies both Wheeler-De Witt 
equation and the boundary condition, namely, 

ipj(a) = Cj<f>j, (35) 

where 4>j is as written above and Cj is again a normalisation constant. 

Note that our wave function is normalised as Jq°° |^| 2 da = 1. Furthermore, we have chosen 
De-Witt boundary (see for example |45|) condition on the wave function. 

iP{a = 0) = 0. (36) 

This simply says that the probability of occurrence of the singularity at a = is zero. We will 
continue to use this boundary condition on the wave function in our following analysis. Existence 
of a normalised wave function with this boundary condition is then suggestive to the fact that the 
universe may not have the initial singularity in the quantum version. Let us also note that the 
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two terms in ip mn individually satisfy the differential equation Hip = 0, but do not individually 
satisfy the boundary condition for M = and Q = 0. With the wave functions at hand, one can 
calculate the mean value of the "radius" of the brane universe, that is 

_ Ja\ip mn \ 2 da 
{a) ~ f\TPmn\ 2 da- (37) 

By dimensional analysis this number can be expected to be the order of We can find the 
solution for T = 0, that is for Ad = — ™. Using the form h(a) as 

2 uwM_ (D - l)u D Q 
a D-2 + _ 2 )a 

Hence, the solution (pj is 

exp[—(j + hiraMp] 



M»)=t-w-s+ r n _x . (38) 



It is important to note that the wave function blows up when the denominator of (|H9|) vanishes for 
non-zero values of a. This happens at the horizon of the bulk black hole. The horizon singularity, 
in black hole physics, can be removed by proper choice of coordinates. Therefore, such divergences 
of the wave function may not be so worrying for us (see |44| for a discussion on this issue). On 
the other hand, the real singularity of the bulk is at a = 0. It is quite easy to check that the wave 
function approaches zero at this singularity satisfying the boundary condition (|36|), 

Let us now proceed to analyse some consequences of the wave function we have just found 
for some specific cases. We first consider the case where M = Q = 0, k = —1 and < 0. The 
classical behavior of the brane was analysed in the previous section. The universe then starts 
from singularity, expands till the maximum size , 1 and then shrinks to singularity. From the 

quantum mechanical perspective, we find the wave function to be 

^ mn ( a ) = ^ _ I e -(m+|>aM p _ e -(n+l)naM p \ (4Q) 



1 + |Ad|o 



4 



Consequently, we can evaluate the expectation value of the universe radius. This turns out to be 

= J jfaq|Vw|^ (41) 

mn | 

Further, expanding the denominator of vpmn and keeping terms to linear in Ad with m = 1, n = 
gives 

1 /ll Ad 395\ m/k ox 

<a> = ^(y + M5Fuij +0(AD) ' (42) 

Similarly, for M^0, Q = 0,k = 1 and Ad = (with T = 0, 1 — >• oo), as written in the second 
line of eq. l(2T]) and is 

^(a) = ^5 _ I e -{m+fraM p _ e -(n+i)naM p \ _ (43) 



4 
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the average value of a for small wbM with m = 1 and n = is 



_jj_ + 3^M r(4-D) /J_ + 1 _ _J_ 
Uw D M r(3 - D) / 1 1 \ 2 



+ 1 -^)+°(^ M ) 2 - ( 44 ) 



We thus notice that in both the above cases, the universe radius gets an average value which, 
in leading order, is proportional to l/M p and is independent of the other dimensionful constants 
that appear in ip m n- One can further evaluate < a > for small but non-zero A. We have explicitly 
checked that the leading order term is again of the order 1/M p and free of other dimensionful 
parameters. This may indicate that the universe, in stead of shrinking to zero size, may stabilize 
at the scale of 1/M p due to quantum effects. In the next section, we would study the solution for 
nonzero T using WKB approximation. 



5 WKB Approximation 



It turns out to be too difficult to find the solution of the Wheeler-De Witt equation for T / 0. 
However, by employing WKB technique, it is possible to get a suggestive picture of the quantum 
universe for branes with non-zero tension. 

Let us recapitulate the recipe to find the wave function of a system using WKB approxima- 
tion. First, we have to solve the the equation H(p, a)ip = Eip, with H(p, a) as the Hamiltonian 
associated to classical system. Setting this classical Hamiltonian, H(p, a) = E, gives us a relation 
between p = f(E, a). Using this approximation, we can find solution to the Hamiltonian equation 
for two regimes: for E > V and E < V, the former corresponds to the classically allowed region 
and the latter to classically forbidden region. The solutions for these cases are 

ipA(a) = = x e ±l f pdx for classically allowed region, 
V("5p) 

ipF(a) = , x e~ -f ^ dx for classically forbidden region. (45) 



Following the prescription as written above, it is quite easy to find the momentum p. It is 
given by 

Ta 



p = M p cosh 1 



= M p ( cosh 1 



(D-l)V^(a) 
Ta 



(D-I)V^HJ 



+ 2mn , 



(46) 



where n is an integer. The conjugate momentum, as written above, has to satisfy the Bohr- 
Summerfield quantization condition 



<j> pda = (m + 8) , 



(47) 
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where m is an integer and associated to "radial" quantum number, and 5 is typically a fraction. 
We will not require the explicit value of 5. The wave function for the classically allowed region is 
then 



Vu(a) 



_L x jj a p d * 



£,—2wnaM p 



x e 



i9(a) 



at 15 " 2 ) ~~ 8(D-2)a^- 4 



A- + a 2 A D + - {D - 1)UJ ° 



where 0(a) is 



9(a) = M p J a cosh' 1 



Ta 



(48) 



(49) 



XD-l)y/K(a). 

It is easy to see from (|48(l that the solution in the classically allowed region is oscillatory in nature 
and simultaneously exponentially damped for n > 0. The solution in the classically forbidden 
region is 

i> F {a) = _L e -2™«M Pe -|e(a)| ; (50) 



where 0(a) is given in eq. ()49|) . 



6 Discussion 



In this paper, we have initiated a systematic study of possible cosmological scenarios on a brane 
where the brane acts as a dynamical boundary of various asymptotically AdS or dS spaces. Our 
emphasize, thought out, was to isolate the cases where big bang singularity was absent on the 
brane. This happens, in particular, when the bulk represents D + 1 dimensional charged AdS or 
dS black hole. Early time singularities for k = 0, ±1 are found to be absent on the brane. This 
was found by explicitly solving the FRW equation when the effective cosmological constant on the 
brane was fine tuned to zero. Furthermore, by studying FRW equation through numerical means, 
we found that the non-singular feature continues even when the effective cosmological constant 
is non-zero. We have also found that, by taking mass and the charge of the bulk black hole to 
large values (but keeping their ratio fixed), the universe can avoid quantum gravity era at early 
time for the models with k = ±1. In section 4 and 5, we analysed the brane universe from the 
quantum mechanical perspective. This was done by employing Wheeler-De Witt equation in the 
mini superspace formalism and also via WKB methods. Within these approximation schemes, we 
studied the wave function of the brane-universe with a special attention to the scenarios where 
classical big bang singularities are present. We then found that the minimum size of the brane is 
of the order of Planck length. 

As it is clear that many interesting issues of cosmological importance were not considered 
in this paper. No matter fields (except those that were induced from bulk) were introduced on 
the brane. These fields would certainly modify the scenarios that we have discussed so far. It 
would be nice to see, however, if the scenario of bounce followed by radiation still persists when 
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such matter fields are included. Furthermore, fluctuation of such matter fields also remain an 
interesting arena to explore. We hope to address these issues in the future. 

Acknowledgment: We wish to thank Rajarshi Ray for discussions. SM would like to thank 
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